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Abstract 

Presence  of  plasma  turbulence  can  strongly  influence  propagation  properties  of  electromagnetic 
signals  used  for  surveillance  and  communication.  In  particular,  we  are  interested  in  the 
generation  of  low  frequency  plasma  turbulence  in  the  form  of  coherent  vortex  structures 
coexisting  with  short  scale  density  irregularities  and  of  lower  hybrid  turbulence.  Lower-hybrid 
type  density  irregularities  are  excited  by  plasma  flows  with  velocity  shear,  whereas  interchange 
or  flute  type  oscillations  in  magnetized  plasma  are  associated  with  Rayleigh-Taylor  type 
instability.  These  types  of  density  irregularities  play  important  role  in  refraction  and  scattering  of 
high  frequency  electromagnetic  signals  propagating  in  the  earth  ionosphere,  inside  a  plasma 
sheath  of  reentry  and  hypersonic  vehicles  and  in  many  other  applications.  We  will  discuss 
generation  of  low  frequency  density  irregularities  due  to  the  presence  of  plasma  flows  with 
velocity  shear  and  interchange  instability. 


1.  Introduction 


The  scattering  of  high-frequency  (HE)  electromagnetic  (EM)  waves  is  a  fundamental 
phenomenon  in  plasma  physics.  In  a  stable  plasma  EM  scattering  occurs  due  to  the  thermal 
electron  density  fluctuations  and  it  is  known  as  “Thomson  scattering”  [1].  Due  to  nonlinear 
interaction  of  waves  in  plasma,  scattering  can  also  arise  from  single  waves  excited  in  the  plasma 
and  is  explained  as  three- wave  nonlinear  interaction  [2],  Scattering  can  also  take  place  due  to 
interaction  of  HF  waves  with  electrostatic  solitons  [3]  or  from  charged  dust  particles  [4].  In  the 
present  paper  we  will  analyze  formation  of  plasma  density  irregularities  due  to  development  of 
lower  hybrid  and  interchange  instabilities.  These  types  of  density  irregularities  play  an  important 
role  in  analysis  of  high  frequency  EM  scattering  and  refraction  in  the  ionospheric  plasma.  The 
reason  for  this  is  that  the  spatial  scales  of  these  plasma  waves  are  comparable  or  smaller  than 
typical  wavelengths  of  EM  signals  used  for  surveillance,  communication  and  OTH  radar 
applications. 
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In  the  terrestrial  high-latitude  (auroral)  magnetosphere  and  ionosphere,  whistler-lower-hybrid 
fluctuations  in  their  electromagnetic  or  electrostatic  limits  are  usually  found  in  association  with 
energetic  electron  beams  along  the  geomagnetic  field  and/or  anisotropic  electron  distributions.  In 
particular,  the  energetic  electron  acceleration  within  the  terrestrial  auroral  region  results  in  the 
emission  of  short-wavelength,  primarily  electrostatic  waves  near  the  resonance  cone  of  the 
whistler-lower-hybrid  branch  [5].  Another  source  of  enhanced  lower  hybrid  waves  in  the  auroral 
ionosphere  is  related  to  strongly  sheared  plasma  flows  driven  by  structured  electric  fields  in  the 
vicinity  of  small  scale  auroral  arcs  [6]. 

At  subauroral  latitudes,  lower  hybrid/fast  magnetosonic  waves  accompany  such  well-known 
space  weather  phenomena  as  subauroral  ion  drifts  (SAID)  and  subauroral  polarization  streams 
(SAPS)  [7,8]  .  Figure  1.1  presents  overview  of  one  of  such  events  occurred  on  18  March  2002 
when  three  independent  spacecraft  (Cluster,  DMSP,  and  Polar)  crossed  the  SAID  channel 
virtually  at  the  same  time  and  within  the  same  magnetic  tube.  It  is  seen  that  the  wave  activity  is 
abruptly  enhanced  inside  the  SAID  channel,  just  interior  to  the  plasmaspere’s  boundary  (the 
plasmapause). 


Figure  1.1  ( left  panel)  World-line  plots  of  the  Cluster,  DMSP,  Polar  satellites  during  the  18-March-2002 
substorm  event.  The  SAID  channels  from  Cluster  1  (Cl)  are  superimposed.  Shown  in  the  middle  are  (top) 
the  C1/C4  outward  electric  field  Ex  in  the  inertial  frame,  (mid)  Cl/EDI  1-keV  electron  counts  (thick  line) 
and  the  cold  plasma  density  no  (thin),  and  (bottom)  C4/STAFF  frequency-time  spectrograms  for  the 
electric  spectral  power  in  (mV/m)2/Hz  (in  log  scale).  Horizontal  lines  indicate  the  lower  hybrid  resonance 
(solid)  and  the  10th,  4th,  and  2nd  harmonics  of  the  H+-ion  gyro  frequency  (dashed),  derived  from  the 
observed  magnetic  field.  Shown  on  the  right  are  the  features  of  the  northern  SAID  channel  from  Polar: 
(top)  the  outward  electric  field,  (mid)  the  plasma  density,  and  (bottom)  spectral  amplitudes  of  electric 
fields  at  32  (dashed  line),  256  (thick),  and  2048  (thin)  Hz. 
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The  remainder  of  this  paper  is  organized  as  follows:  In  Section  2,  excitation  of  the  lower  hybrid 
instability  by  flows  with  velocity  shear  in  the  presence  of  electron-neutral  collisions  is  analyzed. 
Linearized  equations  for  the  lower  hybrid  eigenmodes  inside  a  plasma  slab  with  velocity  shear 
will  be  derived  and  solved  numerically.  Numerical  solution  is  based  on  the  implementation  of 
the  shooting  code.  In  Section  3,  nonlinear  equations  describing  interaction  of  the  lower  hybrid 
waves  inside  a  plasma  flow  with  velocity  shear  with  the  low  frequency  ion  acoustic  oscillations 
will  be  derived,  and  dispersion  equation  for  the  modified  decay  instability  will  be  presented.  In 
Section  4,  the  results  of  the  laboratory  experiments  in  the  NRL  Space  Chamber  on  excitation  of 
lower  hybrid  turbulence  by  ExB  flows  with  velocity  shear  will  be  presented.  In  Section  5,  linear 
and  nonlinear  stages  of  interchange  instability  with  the  ionospheric  plasma  parameters 
corresponding  to  the  equatorial  F-layer  will  be  discussed.  In  the  description  of  the  interchange 
instability  we  retain  finite  ion  Lannor  radius  effects.  Inclusion  of  spatial  scales  comparable  with 
the  ion  Lannor  radius  is  important  for  the  analysis  of  radar  generated  high-frequency  EM  waves 
interacting  with  density  irregularities  associated  with  the  equatorial  spread  F  in  the  ionosphere.  A 
brief  summary  and  applications  are  contained  in  Section  6. 

2.  Excitation  of  Lower  Hybrid  waves  by  a  flow  with  velocity  shear 

In  this  Section  we  will  analyze  instability  of  a  plasma  flow  with  a  transverse  velocity  shear  scale 
length  much  smaller  than  the  ion  gyroradius  but  larger  than  the  electron  gyroradius.  As  shown  in 
[9,10],  under  these  conditions  electrostatic  oscillations  with  the  frequencies  above  the  lower 
hybrid  frequency 


co 


2 

L  H 


1  + 


CO 


CO 


2 

P  e 
2 

c  e 


(2.1) 


can  be  excited.  In  (2.1)  co  and  copj  arc  the  corresponding  electron  and  ion  plasma  frequencies 
and  oo ce  is  the  electron  cyclotron  frequency.  Below  in  the  hydrodynamic  approximation  we  will 
derive  the  equation  for  excitation  of  oscillations  in  the  frequency  range  just  above  the  lower 
hybrid  frequency  coLH .  In  this  equation  we  will  include  wave  dispersion  due  to  the  thermal 
effects,  electron  neutral  collision  frequency  and  nearly  perpendicular  to  the  magnetic  field 

k  ~  TYl 

propagation  wave  vectors  with  —  (m  anc*  M  are  electron  and  ion  masses  respectively),  ky 


-  is  the  wave  vector  perpendicular  to  an  external  magnetic  field  and  kz  -  is  the  wave  vector  of 
oscillations  along  the  magnetic  field.  Electrons  in  the  lower  hybrid  oscillations  are  magnetized, 
but  ions  are  unmagnetized.  The  external  electric  filed  varies  along  the  x-dircction  E0x(x)  and 


creates  electron  flow  with  velocity  shear  directed  along  the  v-direction  V0  (x)  (see  Figure  2.1). 

The  width  of  the  layer  with  electric  field  inhomogeneity  is  smaller  than  the  ion  Lannor  radius.  In 
such  a  system,  we  have  electron  flow  with  velocity  shear  moving  in  y-direction,  whereas  ions, 
being  unmagnetized,  do  not  experience  ExB  drift  and  are  at  rest  in  the  laboratory  coordinate 
system.  We  will  use  linearized  electron  momentum  equation  in  the  drift  approximation  which 
also  contains  electron  neutral  collisions  tenn  and  linearized  ion  momentum  equation  for 
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unmagnetized  ions.  Using  also  continuity  equations  for  electrons  and  ions  and  Poisson  equation 
for  the  electrostatic  potential  of  a  lower  hybrid  wave,  we  can  arrive  to  the  following  linearized 
equation  which  describes  wave  excitation  in  the  non-local  approximation: 

_2  2  _2^  2  2 

3  CO .  ,  ,30  CO  „  CO 

- [(—  (l  +  iTA)rAO  +  ry  - +  v  ((1  + 2-^)1  +  ^)AO  - 

2lx2v  7  z  2  2  7  2  2  y 

Ot  CO  C  Z  CO  co~ 
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(L  +  v  ) 

N  z  en  y 


d2V  30 

- — - 1  +  o)  [L  +v  L  1AO  =  0 

_  9  _  J  pi  L  z  en  z  J 

dx  dy 


(2.2) 


In  equation  (1)  ven  is  an  electron-neutral  collision  frequency,  U0v(x)is  the  flow  velocity  and 

/v  /v  Q  Q 

operator  Lz  is  defined  as:  L _  =  —  +  V0  (x)  —  . 

dt  ’  dy 

We  are  interested  in  the  solution  of  the  equation  (1)  which  represent  waves  propagating  along  the 
y  and  z  directions  inside  a  plasma  slab  and  the  solution  in  the  form  of  the  eigenfunction  in  the  x- 
direction  of  the  flow  inhomogeneity: 

d>  ~  'P(x)  exp(-ico  t  +  k  y  +  k  z)  (2.3) 


Figure  2.1.  In  this  figure  a  nonuniform  electric  field  E{)  (x)  and  an  external  magnetic  field  B(t_  create  a 
nonuniform  electron  flow  with  velocity  shear  along  the  v-direction. 
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We  set  x  —  0  at  the  conducting  surface  of  the  vehicle  and  require 

^(x  =  0)  =  0.  (2.4) 

In  the  free  space  beyond  the  sheath  edge  at  x-L,  the  potential  </>{x)  =  ^(Z)exp[-A:(x-Z)]  so 
that 

=  0.  (2.5) 

Equations  (2.2)  ,(2.4)  and  (2.5)  constitute  the  eigenvalue  problem  to  be  solved  numerically  for 
the  complex  eigenfrequency  co  using  a  standard  shooting  method  [11]. 

Using  equation  (2.2)  for  conditions  used  in  Ganguli  et  al.  and  neglecting  thermal  dispersion  of 
lower  hybrid  waves  (R  =  0),  electron-neutral  collisions  ( ven  =  0 )  and  restricting  analysis  to  the 

waves  with  k_  =  0  ,  the  shooting  code  recovers  the  previously  published 

solutions  co  =  (1.89  +  zl.94)ruii?  [9]From  this  point  using  equation  (2.2),  the  shooting  code  was 

used  to  investigate  the  effect  of  finite  electron-neutral  collisions  on  the  mode.  Figure  2.2  shows 
the  real  frequency  (solid  black  line)  and  growth  rate  (dashed  red  line)  nonnalized  to  the  lower 
hybrid  frequency  as  function  of  electron-neutral  collision  frequency  nonnalized  to  the  electron 
cyclotron  frequency  for  plasma  conditions  similar  to  those  from  published  experimental 
observations  of  this  instability  [12]:  oope/kle  —  8.0,  vE/kleLE  —  13.6,  kyLE  —  0.5,  kzLE  —  0, 
and  argon  ions.  The  following  electric  field  profile  was  used  in  the  code: 


E0x(x)  =  E0  sec  h2 (-^-)  (2.6) 


where  LE  is  the  velocity  shear  scale  length  in  the  electron  flow  and  the  shear  flow  velocity  along 
the  y-direction  is  defined  as: 


V0v(X)e  =  c 


E0(x)xB0 

Bl 


(2.7) 


where  c  -  is  the  speed  of  light. 

Results  of  numerical  solution  of  equation  (2.2)  with  inclusion  of  electron-neutral  collisions,  but 
with  omission  of  thermal  dispersion  and  with  kv  =  0  are  presented  by  Figure  2.2.  It  can  be  seen 
that  presence  of  weak  collisions  has  a  strong  effect  on  both  the  growth  rate  and  real  frequency, 
though  the  mode  is  not  completely  stabilized.  We  stop  at  ven~0.1o)ce  to  remain  in  the  regime 
where  the  drift  approximation  is  valid.  However,  it  is  expected  that  ven~D.e  is  needed  before 
collisions  could  totally  suppress  excitation  of  the  mode. 
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Figure  2.2:  Normalized  real  frequency  (solid  black  line)  and  growth  rate  (dashed  red  line)  dependence  on 
normalized  collision  frequency. 

To  obtain  solution  for  the  frequency  and  the  growth  rate  of  excited  modes  in  the  local 
approximation,  we  again  will  use  equation  (2.2)  neglecting  thermal  dispersion  of  lower  hybrid 
waves(R=0),  electron-neutral  collisions  ( ven  =  0  )  and  restricting  analysis  to  the  waves  with 

k_  =  0  .  This  will  allow  us  to  examine  the  qualitative  behavior  of  this  instability  without 
employing  the  full  numerical  solution.  To  arrive  at  the  local  approximation  for  the  lower  hybrid 
modes,  we  start  from: 


where 


—  k2  <p  —  0, 


_ 2  -2  8 2  m2  kyvE 

V  +  1  hJ2  —  1  th  —  kyVE 


and  8  —  cope/{2e,  ky  —  kyLE,  go  =  gj/gl>lh,  and  vE  —  ve/g)lhLe.  We  obtain  the  local 

dispersion  relation  by  taking  k2  —  0  and  using  the  values  for  vE  and  vE  at  x  —  0  as  in  Romero  et 
al.  [9].  This  results  in  a  cubic  equation  for  m: 

_3  (  2  82  \_2  _  __ 
a)  +  kyvE o  \j^S2  -u  +  kyvEO  =  0, 

Figure  1.3  shows  the  local  and  non-local  solutions  for  the  following  parameters:  8  —  8.0, 
vE  —  1.1,  kyLE  —  0.5,  kzLE  —  0,  and  argon  ions.  The  local  solution  does  a  good  job  of 
approximating  the  general  behavior  of  the  non-local  solutions,  but  more  importantly,  it 
represents  a  good  operating  regime  for  studying  the  non-linear  behavior  of  this  mode.  The 
maximum  growth  rate  occurs  for  a  real  frequency  just  above  the  lower  hybrid  frequency,  and  the 
mode  only  exists  for  a  small  band  of  perpendicular  wave  numbers. 
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Figure  1.3:  Normalized  real  frequency  (solid  black  line)  and  growth  rate  (dashed  red  line)  dependence  on 
normalized  perpendicular  wave  number  comparing  local  and  non-local  solutions. 


3.  Nonlinear  interaction  of  lower  hybrid  and  ion  acoustic  waves 


It  is  well  known  that  the  most  effective  mechanism  for  saturation  of  the  lower  hybrid  instability 
is  connected  with  the  decay  and  modified  decay  instabilities  which  occur  due  to  interaction  of 
the  excited  lower  hybrid  waves  with  the  ion  acoustic  oscillations  [13,  14].  In  the  experiments 
described  in  Section  4,  the  electron  temperature  was  much  larger  than  the  ion  temperature  with 
T 

~  10 .  This  justifies  implementation  of  the  nonlinear  scheme  with  involvement  of  ion  acoustic 

waves  instead  of  the  nonlinear  interaction  scheme  based  on  induced  scattering  on  electrons. 
System  of  equations  for  the  electrostatic  potential  of  lower  hybrid  waves  excited  in  a  system 
consisting  of  unmagnetized  ions  and  moving  through  them  in  y-direction  magnetized  electrons 
with  velocity  shear  in  x-direction  with  the  speed  V0y  (x)  and  nonlinearly  coupled  with  ion 

acoustic  oscillations  is  presented  below.  Density  perturbations  in  ion  acoustic  oscillations  are 
defined  as  Sn  .  This  system  of  equations  has  the  following  form: 


92  co  .  9  't  o  a2® 

— [(^(l  +  f?2A)T;A®  +  <  — 
9t  co„  9  z 


col  7  d2V0y  a® 


CO 


f L. 


9  x2  9  y 


]  + 


9  ~9  CO„„ 

C o2piL2zm  =  -f 


cot  9t ' 


■LA— .«}  (3-D 

na 


^-(Kr;  +  r>  = 

9t2  AnM  co . oo, 


-  A{®,®*} 


(3.2) 


8ci  db  8a  8b 

In  equations  (3.1)  and  (3.2)  Poisson  brackets  {a,b}  are  defined  as  {a,b}  = - and 

9x  9y  9y  9x 

L  • 

z  dt  °>A  ay 

Below,  in  the  local  approximation  we  will  derive  dispersion  equation  for  modified  decay  of 
excited  lower  hybrid  waves  and  find  the  estimate  for  the  threshold  of  this  instability.  In  the 
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absence  of  the  electron  flow,  modified  decay  of  magnetosonic  waves  in  the  vicinity  of  the  lower 
hybrid  frequency  was  analyzed  in  [15]. 

We  will  examine  the  case  when  the  pump  wave  is  excited  due  to  the  presence  of  velocity  shear 
in  the  electron  flow.  The  parametric  decay  of  a  pump  wave  into  two  satellites  in  vicinity  of  the 
lower  hybrid  frequency  and  low-frequency  ion  acoustic  type  mode  represented  by  the  low- 
frequency  density  fluctuations  can  be  describe  as  follows.  We  will  introduce  the  pump  lower 
hybrid  wave: 


'  /«k0  \  i(kor-«V) 


k°)e 


+  c.c  =  '¥u  /k°r"v) 


+  c.c 


As  well  as  Anti-Stokes: 


1  iak[  .  i(k,r-a^1t) 


®k,  =  (  ^  ‘Pk,  e  ^  )e 


+  c.c  =  '¥k  eHk'r-c%,) 


+  C.C 


and  Stokes  satellites: 


(3.3) 


(3.4) 


1  ‘<*k2  \  i(k2r-«i20  .  „„  IT/  J(k2r-®k2<) 


®k2  =(<2(p^e  k2)e 


» T(  L\ 

+  c.c  =  ^k  e 


+  c.c 


(3.5) 


For  the  ion  acoustic  type  density  fluctuations  we  will  use: 


1 


SnK  =  {-8nke^)e 


iak  n  /(kr+K  i)  r  /(kr+vk  t) 

1  +c.c  =  SN,e  1 


k» 


+  C.C 


(3.6) 


The  following  conditions  on  the  wave  vectors  of  two  satellites  kj  and  k2  as  well  as  the  wave 
vector  in  ion  acoustic  perturbations  ksare  imposed:  kl  =k0  +ks  and  k2  =k0  -ks.  Now,  using 
equations  (3.1)  and  (3.2)  the  following  dispersion  equation  for  the  modified  decay  instability  of 
the  lower  hybrid  waves  interacting  with  low  frequency  ion  acoustic  density  fluctuations  can  be 
derived: 


2  2 

•2-r^2  'T^'-21  ^ - — ^2[koxkJ2— ®2|Tkof  X 

47TMn()  (DW„  (Q  „ 


V  =WT)  +  v;)k;}~ 


,  1  [(^®k0  +  ksy  Km  +  V  )i  1  (<K0  “  ksy  Km  ~V), 

x{ — - - 2 - - - + - 2 - - - } 

D«k,  i(®k  -®k2)  +  »']  DRk,  i(®k  -®tl)-V] 


To  obtain  (3.7),  we  also  used: 

=  ^lh  +  ,  <x>LH  ~  k0yV0m, 


dXy  2JA 


dx1 


L- 


(3.7) 
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Finally,  from  equation  (3.7)  the  following  estimate  for  the  threshold  of  the  modified  decay 
instability  of  the  lower  hybrid  waves  can  be  obtained: 


(3.8 

4/r«07;  Mco2pe 

Using  (3.8)  it  is  also  possible  to  estimate  the  amplitude  of  the  low  frequency  ion  acoustic  type 
density  fluctuations  which  appear  in  the  system. 

4.  Experiments  on  excitation  of  lower  hybrid  turbulence 
in  the  NRL  Space  Chamber 

In  a  series  of  experiments,  Amatucci  et  al.  [12]  demonstrated  that  strongly  sheared  electron  flows 
perpendicular  to  the  background  magnetic  field  with  scale  size  smaller  than  an  ion  gyroradius  but 
larger  than  an  electron  gyroradius  can  drive  electrostatic  oscillations  near  the  lower  hybrid 
frequency.  These  observations  were  consistent  with  the 

theoretical  predictions  of  Ganguli  et  al.  [8]  and  Ganguli  and  Romero  [9]  and  lent  s  upport  to  the 
relaxation  scenario  for  highly  stressed  magnetospheric  boundary  layers.  For  example,  during 
intense  solar  activity,  the  plasma  sheet  boundary  layer  can  become  highly  compressed,  with 
gradients  in  the  plasma  across  the  layer  self-consistently  generating  localized  electric  fields  and 
highly  sheared  flows  across  the  magnetic  field.  These  sheared  flows  generate  plasma 
instabilities  that  work  to  dissipate  differential  flows  and  widen  the  boundary  layer.  An 
analogous  situation  can  arise  in  the  flowing  plasma  surrounding  a  hypersonic  vehicle,  resulting 
in  the  generation  of  strong  plasma  turbulence. 

The  experiments  on  the  generation  of  shear-driven  lower  hybrid  waves  were  conducted  in  the 
Space  Physics  Simulation  Chamber  (SPSC)  at  the  Naval  Research  Laboratory.  The  main 
chamber  part  of  the  SPSC  consists  of  a  1.8  m  diameter,  5  m  long  cylindrical  vacuum  vessel 
designed  for  controlled,  scaled  experiments  of  space  plasma  processes.  For  these  experiments,  a 
large  area  (~  75  cm  diameter)  hot  filament  plasma  source  was  placed  at  one  end  of  the  chamber, 
while  a  smaller  filament  source  (~  2  cm  diameter)  is  placed  at  the  opposite  end.  As  can  be  seen 
in  a  schematic  of  the  experimental  setup  in  Figure  4.2,  an  isolated  electrode  blocks  a  small 
portion  of  the  plasma  generated  by  the  large  area  source,  which  is  filled  with  the  plasma 
generated  by  the  small  plasma  source.  The  blocking  disk  and  the  outer  grid  of  the  small  source 
allow  for  control  over  the  plasma  potential  of  the  small  area  plasma  column.  Strong  radial  dc 
electric  fields  can  be  generated  by  this  steep  variation  of  plasma  potential  across  the  cylindrical 
boundary  layer  between  the  two  plasmas.  This  in  conjunction  with  the  filament  heater  currents 
on  both  sources  enables  independent  control  over  the  electric  field  and  density  gradients  across 
the  boundary  of  the  two  interpenetrating  plasmas,  simulating  the  stressed  plasma  sheet  boundary 
layer  conditions. 
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Figure  4.2:  Schematic  diagram  of  the  experimental  setup  showing  the  two  interpenetrating  plasmas. 

Typical  parameters  for  the  steady-state  argon  plasma  are:  plasma  density  n  ~109-1010  cm'3,  ion 
and  electron  temperatures  7}  -0.05  eV  and  Te  -0.25  eV,  ion  and  electron  thermal  speeds  vti  - 
5X104  cm/s  and  vte  ~  2x1 07  cm/s,  and  an  axial  magnetic  field  B  =  40  G.  This  yields  ion  and 
electron  gyrofrequencies  £2,  -1.5  kHz  and  Qe  -110  Mhz,  ion  and  electron  gyroradii  p,  -3.5  cm 
and  pe  -0.03  cm,  and  lower  hybrid  frequency////  -400  kHz.  The  neutral  density  nn  is  variable 
from  1 0 1 3- 1 0 14  cm'3,  and  the  plasma  column  diameter  and  effective  length  are  75  cm  and  3  m, 
respectively. 

Plasma  potential  was  measured  using  radial  emissive  probes.  The  derivative  of  the  resulting 
profile  is  the  electric  field  profile.  The  measured  electric  field  yields  a  typical  transverse  electric 
field  scale  length  LE  from  0.6- 1.0  cm  (0.17-0.3  pt)  and  the  magnitude  can  be  controlled  up  to  40 
V/cm.  Typical  plasma  potential  (solid  line  without  symbols)  and  electric  field  (dashed  line) 
radial  profiles  can  be  seen  in  Figure  4.3.  Since  the  electric  field  scale  length  is  in  the  range  pe  < 
Le  <  pi,  the  electrons  ExB  drift  but  the  ions  do  not.  For  sufficiently  strong  electric  fields,  an 
instability  with  frequency  in  the  lower  hybrid  frequency  range  is  observed  within  the  shear  layer. 
Figure  4.3  shows  an  overlay  of  mode  amplitude  (solid  circles)  as  a  function  of  radial  position 
localized  to  the  edge  of  the  transverse  electron  flow  layer.  A  sample  wave  spectrum  is  also  inset 
showing  a  mode  with  frequency  -780  kHz  (-1.9 fui),  indicating  that  sheared  transverse  electron 
flows  without  the  presence  of  a  density  gradient  can  drive  lower  hybrid  waves. 


10 

Approved  for  public  release;  distribution  unlimited. 


radial  position  (cm) 

Figure  4.3:  Radial  profile  of  lower  hybrid  wave 
amplitude  for  the  uniform  density  case  and  wave 
spectrum  as  inset. 


Figure  4.4:  Comparison  of  experimentally  observed 
(filled  symbols)  and  theoretically  predicted  (line) 
values  of  mode  frequency  as  a  function  of  magnetic 
field  strength. 


Figure  4.4  shows  a  comparison  between  theoretically  predicted  mode  frequency  for  these 
velocity  shear-driven  waves  and  the  experimentally  measure  values  as  a  function  of  the  magnetic 
field  strength.  There  is  good  agreement  between  the  theoretical  predictions  and  the  measure 
mode  frequencies.  The  observed  instability  also  exhibits  a  threshold  electric  field  (~7.5  V/cm) 
for  the  waves  to  appear.  This  coupled  with  the  spatial  localization  of  the  mode  amplitude  to  the 
shear  layer  suggests  that  the  shear  in  the  electron  ExB  flow  is  responsible  for  driving  the 
observed  lower  hybrid  waves. 


While  these  experiments  were  scaled  to  magnetospheric  boundary  layer  conditions,  they  could 
equally  be  applied  to  hypersonic  vehicle  plasmas.  Future  experimentation  will  focus  on 
increasing  the  overall  size  of  the  flowing  plasma  region  while  maintaining  the  steep  gradients. 
By  driving  the  mode  well  in  the  nonlinearly  saturate  state,  we  will  investigate  the  formation  of 
turbulent  plasma  density  structures  and  the  effect  of  such  structures  on  the  scattering  and 
transmission  of  electromagnetic  waves. 


5.  Interchange  instability  in  the  ionospheric  F-layer 

Below  we  will  present  analysis  of  generation  of  density  vorticies  on  the  nonlinear  stage  of 
interchange  instability[16  -18].  Such  density  vorticies  can  play  significant  role  in  scattering  of 
high-frequency  EM  waves  generated  on  the  nonlinear  stage  of  excited  instability.  Such  vorticies 
appear  in  the  ionospheric  F-layer  and  are  part  of  the  ionospheric  clutter  [19  -  22],  They  also  can 
exist  in  a  plasma  flow  around  a  hypersonic  vehicle  and  adversely  effect  sensors  performance. 
This  study  is  also  relevant  to  a  laser  probing  in  a  high  energy  density  plasmas  and  for  diagnostics 
of  the  edge  plasma  in  tokamaks.  Our  analysis  will  be  applicable  to  both  low-beta  and  high-beta 
plasmas  and  will  include  finite  ion  Lannor  radius  effects  [16,  17].  Inclusion  of  spatial  scales 
comparable  with  the  ion  Larmor  radius  is  important  for  the  analysis  of  radar  generated  high- 
frequency  EM  waves  interacting  with  density  irregularities  associated  with  the  equatorial  spread 
F  in  the  ionosphere. 
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5.1  Basic  Equations 


Two-fluid  macroscopic  equations  will  be  used  to  describe  low-frequency  interchange  modes  ( 
where  a>  is  the  frequency  of  the  flute  mode  and  Q.  =  ^  is  the  cyclotron 


frequency  of  the  ion  with  charge  Z  and  mass  mi  =  /urn p ,  with  mp  being  the  proton  mass  and  c 
is  the  speed  of  light)  in  a  weakly  inhomogeneous  plasma  with  external  magnetic  field  B0,(x)  ■ 
As  is  customary  in  the  interchange  turbulence,  oscillations  are  taken  to  be  unifonn  in  the 
direction  of  the  magnetic  field,  i.e.,  the  wave  vector  along  the  magnetic  field  k}l  =  0 .  We 

consider  a  weakly  inhomogeneous  high-beta  plasma  of  slab  geometry  with  density  n0(x)  in  the 
presence  of  an  inhomogeneous  magnetic  field  B0z(x) ,  where  x  is  the  direction  of 
inhomogeneity.  We  will  also  retain  the  gravity  term  g  =  ger ,  which  drives  the  instability  of  the 

flute  modes.  One  of  the  main  difference  if  to  compare  with  the  case  of  a  low-beta  plasma  is 
that  in  a  high-beta  plasma  the  electric  field  in  the  flute  oscillations  is  not  irrotational  (i.e. 

VxE^O)  as  in  the  low-beta  case  and  is  written  as  E  = -V  x  ®  -  ( J/)(dA/^) .  The  plasma 

density  N  can  be  expressed  as  the  sum  of  a  slowly  varying  component  with  x,  the  equilibrium 
plasma  density  n0(x)  ,  and  a  perturbed  component  due  to  interchange  oscillations  5n(x,y,t )  . 
Likewise,  the  magnetic  field  is  written  as  B,  =  B()  ,  +  8 B  (x,y,t) .  We  also  assume  the  quasi¬ 
neutrality  condition  Ne  =  ZNf .  The  equilibrium  condition  for  the  case  under  consideration  can 
be  written  as 


where 


P  „  .  gfi, 


KR=—K„+-  1 

B  2  "  2  Vi 


(5.1) 


1  dn0 
n0  dx 


>0, 


kh 


B0:  dx  50- 


Parameter  /?  -  ion  plasma  beta,  defines  the  ratio  of  the  kinetic  energy  containing  in  plasma 


ions  to  the  magnetic  field  energy.  Below  we  will  also  use  parameter  /?,  = 


Bl 


-  electron 


plasma  beta  and  total  plasma  beta  /?  =  /?.+  (de . 

To  obtain  the  coupled  nonlinear  system  of  equations  for  density  N,  electrostatic  potential  <f>, 
and  magnetic  field  B_ ,  we  start  from  the  quasi-neutrality  condition 


<i/vj  =  0  (5.2) 

The  continuity  equation  for  the  electron  density 
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-  +  div(N  v  )  =  0 


And  the  equation  for  the  magnetic  field  in  the  flute-type  oscillations 

^  =  —  h  (5-4) 

dt  c 

To  the  system  (5.2)  -  (5.4)  we  should  add  equations  of  motion  for  ions  and  electrons.  The 
equations  of  motion  for  the  ions  are  written  in  the  fonn 

+  -V)v,]  =  qiNj(E  +  — [v,  x  (B0  +  AB)])  -  XP.  -X  ft  +  g  (5.5) 

ot  c 

where  qt  =  Ze ,  and  Pt  =  NT .  It  is  also  assumed  that  the  viscosity  tensor  ft  contains  only 
gyroviscosity  components 


N  T  dv  dv 

n  =  -n  = - — ( — —  h - L) 

""  2Q,  dy  dx 


N,T,  dvx 

n  =n  = — — -( — - — ) 

*v  “  2Q,.  dx  dy  (5  y) 

Below  it  is  assumed  that  the  electric  drift,  the  diamagnetic  drift,  and  the  gravitational  drift  are 
all  of  the  same  order,  and  of  first  order  in  the  small  parameter  s  defined  by: 

_  Zf®.  _  51jl  Sni  K'«  kb  fkp  f  e 
coci  T  B0:  n0i  k  k  '  (58) 

From  the  ion  equation  of  motion  one  can  obtain  an  iterative  solution  for  the  ion  velocity: 

c  _  vl  v«..  p  i  r  dv!n)  v-n.w 


V/"+1)  =  -  — zxE  +  -^zx^--zx  —  zx^—  +  zx(v,(")-v)v,("] 

B0:  coci  /I,  ojci  coci  dt  1  ' 


Next,  we  can  estimate  the  order  of  various  terms: 


m'n'0J“  (5.9) 


c  1  SA  c  coA 

VF  =  —  zxV®  +— zx - Mh - 

B  B  dt  B  B 


V„,  = — -zx 


c  coA  ,  J Ze®]  VT  ( 8B](  a>  \  , 

m+  VTi  (kPi)  —  +-J-  —  —  ~VTi{kPi)i 
B  B  \  Tj  )  kpi  v  B 

(  K"  f)  yi  ^ 

V  In  ni  ~  VTi  {kPt )  -f  +  ^  ~VTi{kPi)s 

\  k  n0i  ) 


X  g  =  -zx 


g  8 


°>  a  ®  a 


V T i  (Ptk)-£~~  V T i  (Pik)S 


(5.10) 

5.(11) 

(5.12) 
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The  zero  order  velocity  is  the  equilibrium  diamagnetic  drift  plus  the  gravitational  drift 
(assuming  no  equilibrium  electric  field,  E0): 


y_(°)  =  2kzxZ"ro_zxJL  =  v  +v 

00,, 


(5.13) 


The  first  iteration  gives: 


d1)  _  c 


- z  x 

Bn 


V  nj 

x - - 


v  n 


<0) 


Z  X  ■ 


—  Z  X  ■ 


CO,;  11; 


CO „ 


(5.14) 


=  V  +  V  +  V  +  V 

’  e  T  T  Di  T  ’g  *no> 

do) 


where  V-I1<0)  is  computed  at  V)'>! .  And  from  the  next  iteration  one  obtains  the  following 
solution: 


\42)  =-— zx£E+-^-zx^-zx-^-+— 


Bn 


(»a  «,  0>ci  <°ci 


ZX - - - hzx 

dt 


•  v)v41] 


v-n(1) 

+zx-— —  (5.15) 


where  V  -IIll)  is  computed  at  V. 


(i) 


5.2  Corrections  to  Gyroviscous  Cancellation 

In  the  case  of  low-beta  plasma  the  gyroviscous  flux  cancels  the  diamagnetic  flux  (see  for 
example  [23]  and  [24])  .  Below  we  will  show  that  in  a  high-beta  plasma  this  cancellation  does 
not  take  place  and  as  a  result  new  terms  appear  in  the  final  nonlinear  equations  that  describe 
interchange  turbulence.  Using  equation  (5.5)  we  can  write  the  expression  for  the  perpendicular 
ion  velocity  keeping  only  the  terms  that  are  first  order  in  s : 

Vu=Vs+V„+VJ+0(«2)  (5.16) 

For  the  electric  drift  VE ,  the  diamagnetic  drift  VWi  and  for  the  gravitational  drift  Vg  we  have  the 
following  expressions: 


V£  =--C-zx5E  VDi  =  — 

B0  ooci  n,  coci 


Now  the  ion  continuity  equation  can  be  written  as: 

"5V, 


-» 

dt  co  - 


V< 


dt 


+(v*-v)vL 


11  w 

f  z  0V„.  I 

+v- 

X  Dl 

JI 

0e 

CD 

+ 


— +(v„.v)vj] 

[COci 

L  mini  JJ 

=  0 


(5.17) 


(5.18) 


14 

Approved  for  public  release;  distribution  unlimited. 


After  lengthy  calculations  we  can  arrive  at: 


a>ci  mtnt  2  a)ciB0z  Gt  coci  dy 


(5.19) 


In  a  low-beta  plasma  relation  (5.19)  should  be  replaced  by  the  well  known  gyroviscous 
cancellation  relation  (see  for  example  [23]  and  [24]): 


V-Ax[— +  (V„-V)VJ}  =  0 


ao,,  mini 


(5.20) 


5.3  Nonlinear  Equations  for  Description  of  Interchange  Turbulence 


Finally,  with  the  help  of  relation  (5.19)  we  can  arrive  at  the  following  complicated  system  of 
nonlinear  equations  for  density  N,  electrostatic  potential  ©  and  magnetic  field  SB,  : 


2  ,4 

dA±Q  V-  kn  +  kb  GA  ©  |  gB0i  lPl  2  A  , dSn 
dt  coci  l  +  dy  ni0c  1  +  4  1  Gy 

2  2 


VG  (8A±SB:  |  g  dAJB  g  {dA±0  |  V*  B0:  8 A  Jn) 
cocic  dt  mci  dy  mci  dy  cocic  ni0  dy 


(5.21) 


Vl 


o)cin0i 


-V±-{V±0,Sn}  +  —  {  A±©,©} 


B, 


0z 


dSn 

dt 


c  kn+kb  5© 

+  K  i+4  ty 
2 


PXlKJ_  2KN8Sn 

Pi  0)n  i+4  fy 

2 


5B_ 


(5.22) 


(5.23) 


In  equations  (5.21)  -  (5.23)  we  used  quasi-neutrality  condition  Sni=Sne=Sn  and  Poisson 
brackets  are  defined  as: 


{a,b} 


da  8b  da  db 
dx  dy  dy  dx 


Equations  (5.21)  -  (5.23)  describe  the  excitation  and  nonlinear  evolution  if  compressible 
electromagnetic  interchange  modes  in  a  finite  beta  plasma  with  inhomogeneous  density  and 
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magnetic  field.  As  one  can  see,  the  nonlinear  terms  in  equations  (5.21)  -  (5.22)  are  rather 
complicated.  They  contain  so-called  vector  nonlinearities  (represented  through  Poisson 
brackets)  which  are  the  source  for  generation  of  large  scale  vortex  structures  coexisting  with 
short  scale  spectral  components  produced  on  the  nonlinear  stage  of  modulation  instability  of 
interchange  modes.  Equations  (5.21)  -  (5.23)  will  be  used  to  analyze  interchange  turbulence 
and  scattering  of  high  frequency  electromagnetic  waves  on  density  irregularities  associated 
with  interchange  turbulence. 

5.5  Instability  of  Interchange  Modes 

The  dispersion  equation  for  interchange  mode  instability  can  be  obtained  from  the  linearized 
system  of  equations  (5.21)  -  (5.23).  The  dispersion  relation  was  solved  numerically  for  the 
following  set  of  plasma  parameters,  roughly  approximating  the  plasma  conditions  that  exist  in 
the  so-called  spread-F  layer  in  the  ionosphere: 

n0  ~  2  x  103  cm  3  -  plasma  density 

Bu  ~  0.4  G  -  magnetic  field 

g  ~  500  cm  Is 2  -  gravitational  constant 

Tt  =  Te  ~  1.4  eV  -  ion  and  electron  temperature 

mi  =  16  x  mp  -  ion  mass  and  mp  -  proton  mass 

Results  of  the  numerical  solution  of  the  linearized  dispersion  equation  are  presented  in 

Fig.  5.1  below.  In  the  limit  of  very  small  [3  ,  the  linear  dispersion  relation  can  be  solved  and  the 

linear  growth  rate  should  scale  as: 

/max=V^g-  (5-24) 

In  Fig.  5.1  we  compare  the  growth  of  the  energy  in  the  density  and  potential  fluctuations  from 
the  numerical  solution  of  equations  (5.21  -5.23)  with  the  maximum  growth  rate  calculated  from 
Eq.  (5.24).  From  Eq.  (5.24),  we  find  that  ymax  ~  0.0408  A1 ,  or  ymax  /  coci  ~  1.7 xlO’4 .  As  shown  in 

Fig.  5.1,  we  find  that  ymax  /  coci  ~  1.6x10  4  provides  a  reasonable  fit  to  the  code  results. 


16 

Approved  for  public  release;  distribution  unlimited. 


IE-5 


Figure  5.1.  Density  and  potential  growth 


5.6  Numerical  Solution  of  Nonlinear  Equations 

Below  we  present  results  of  numerical  solution  of  nonlinear  system  of  equations  (5.21)  -  (5.23). 
For  the  numerical  analysis  we  use  the  code  Flute  developed  by  Dr.  J.N.  Leboeuf;  detailed 
description  of  the  code  can  be  found  in  [17]. 

The  linear  growth  of  the  electrostatic  wave  energy  and  density  perturbations  followed  by  the 
nonlinear  saturation  of  the  instability  is  presented  in  Fig.  5.2. 


Figure  5.2.  Linear  growth  and  nonlinear  saturation  of  the  interchange  instability 
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Snapshots  of  the  density  fluctuations  are  shown  Fig.  5.3.  The  initial  perturbations  applied  to  the 
vorticity  and  density  are  random,  with  the  normalized  amplitude  10-4.  This  is  illustrated  in  Fig. 
5.3(a).  The  linear  phase  of  the  instability  is  illustrated  in  Fig.  5.3(b)  and  5.3(c).  It  is  clearly  see 
that  perturbations  are  elongated  along  the  x  -  direction.  This  are  so-called  streamer-like 
perturbations.  The  transition  from  linear  to  nonlinear  and  the  nonlinear  phase  is  captured  in  Fig. 
5.3(d)  -  (f).  On  the  nonlinear  stage,  when  amplitudes  of  perturbations  exceed  the  threshold  for 
the  modulation  instability,  transition  from  streamer-type  perturbations  to  the  zonal  flow  type  of 
perturbations  (elongated  along  the  y  -  direction,  perpendicular  to  the  direction  of  plasma 
inhomogeneity)  takes  place.  This  is  one  of  the  expected  characteristics  of  the  interchange 
turbulence.  The  nonlinear  phase  is  also  characterized  by  a  broad  spectrum  of  excited 
wavelengths  and  frequencies. 


(a)  Qj  =  0 


( b )  n(f  =  10 


(c)  Q,t  =  300 


Density  in  linear  stage 


(d)  =  2700 


(e)  Q,i  =  2800 


(f)  Cltt  =  2980 


Density  in  nonlinear  stage  and  formation  of  zonal  flows  and  nonlinear  vortex  structures 

Figure  5.3.  Snapshots  of  the  density  fluctuations 
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6.  Summary 


This  paper  consists  of  two  parts.  In  the  first  part  (Sections  2-4)  results  of  the  theoretical 
analysis  and  the  NRL  Space  Chamber  experimental  study  of  excitation  and  nonlinear  saturation 
of  lower  hybrid  instability  in  plasmas  with  unmagnetized  ions  and  magnetized  electrons  drifting 
across  the  magnetic  field  due  to  the  presence  of  an  external  inhomogeneous 
electric  field  E0x(x)  are  presented.  Drifting  electron  flow  has  velocity  shear  and  is  the  source  of 

the  instability  (see  [8,9]  for  details).  Equation  for  excitation  of  instability  in  the  vicinity  of  the 
lower  hybrid  frequency  with  inclusion  of  the  thermal  dispersion  and  electron-neutral  collisions 
was  derived  and  solved  in  nonlocal  and  local  approximations.  Results  show  satisfactory 
agreement  for  both  numerical  solution  obtained  with  the  help  of  the  shooting  code  and  solution 
in  the  local  approximation.  Next,  nonlinear  equations  describing  coupling  of  excited  lower 
hybrid  waves  with  low  frequency  ion  acoustic  perturbations  were.  Using  these  equations, 
modified  decay  instability  of  the  lower  hybrid  waves  in  a  local  approximation  was  analyzed.  In 
particular,  dispersion  equation  for  the  modified  decay  instability  was  derived  and  its  threshold 
was  determined. 

In  the  second  part,  analysis  of  interchange  instability  with  inclusion  of  finite  ion  Larmor  radius 
effects  in  a  high  beta  plasma  is  presented.  Obtained  results  are  used  to  analyze  linear  and 
nonlinear  stages  of  the  instability  as  well  as  formation  of  vortex  structures  in  a  plasma  with 
parameters  similar  to  that  in  the  ionospheric  F  layer.  Numerical  results  show  that  formation  of 
vortex  structures  occurs  when  the  amplitude  of  density  perturbations  exceeds  the  threshold  for 
the  modulation  instability  and  on  the  nonlinear  stage  density  perturbations  in  the  vortex 
structures  can  be  very  large.  Developed  technique  can  be  also  applied  to  analyze  interchange 
modes  in  high  energy  density  physics,  such  as  Z-pinch  plasmas.  This  study  is  also  relevant  to  a 
laser  probing  in  a  high  energy  density  plasmas  and  for  diagnostics  of  the  edge  plasma  in 
tokamaks.  Another  areas  of  applications  -  magnetospheric  and  space  plasmas. 

The  types  of  density  irregularities  examined  in  the  paper  play  important  role  in  refraction  and 
scattering  of  high  frequency  electromagnetic  signals  propagating  in  the  earth  ionosphere,  inside  a 
plasma  sheath  of  reentry  and  hypersonic  vehicles  and  in  many  other  applications  related  to 
surveillance  and  navigation.  Detailed  analysis  of  high  frequency  electromagnetic  wave  refraction 
and  scattering  in  the  presence  of  developed  low  frequency  plasma  turbulence  described  in  the 
paper  is  the  subject  for  the  future  work. 
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